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We introduce a generalized method of holonomic quantum computation (HQC) based on encoding 
in subsystems. As an application, we propose a scheme for applying holonomic gates to unencoded 
qubits by the use of a noisy ancillary qubit. This scheme does not require initialization in a subspace 
since all dynamical effects factor out as a transformation on the ancilla. We use this approach to 
show how fault-tolerant HQC can be realized via 2-local Hamiltonians with perturbative gadgets. 
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Introduction. — A basic requirement for the construc- 
tion of a reliable quantum computer is the ability for 
high-fidelity storage and manipulation of quantum infor- 
mation. Storage is most generally achieved through en- 
coding in suitable physical degrees of freedom that may 
be protected from unwanted interactions through sym- 
metries or by active operations. At the same time, ma- 
nipulating quantum information requires addressing the 
relevant degrees of freedom in a precise and robust way. 

Holonomic quantum computation (HQC) is one 
method that promises a resilient way of information pro- 
cessing through all-geometric, adiabatic control [1]. In 
this approach, logical states are encoded in a degener- 
ate eigenspace of the Hamiltonian, and gates are realized 
by adiabatically varying the Hamiltonian along suitable 
paths in the space of control parameters. This gives rise 
to geometric transformations inside the eigenspace Q, 
which depend on certain global properties of the path 
and are thus robust against local fluctuations that pre- 
serve those properties [3| . This method employs encoding 
in a subspace, which is only a special case of the most 
general form of encoding possible [3] — encoding in virtual 
subsystems [11 . The latter type of encoding has numerous 
applications in the area of decoherence control, ranging 
from noiseless subsystems that offer protection through 
more efficient encoding or in cases where no subspace 
protection exists 0], to operator error-correcting codes 
[7| that allow for simplified recovery methods leading to 
improved fault-tolerance thresholds 3. Given the oper- 
ational robustness of the holonomic approach, a natural 
question is whether a generalized method for HQC com- 
patible with encoding in subsystems is possible. 

In this paper, we answer the above question affirma- 
tively. The paper reports three main results. The first 
one is a general framework for HQC in subsystems. A 
distinctive feature of this framework is that it involves 
performing simultaneous computations in the different 
eigenspaces of the Hamiltonian. We show that given suf- 
ficient control over the parameters of a Hamiltonian, it is 
possible to generate any combination of geometric trans- 
formations in its eigenspaces. This possibility has been 
suggested in Refs. [5|, |9|, but a proof has been lacking. 
A remarkable consequence of this result is the possibility 
to apply purely geometric operations on a given system 
by the use of a noisy ancilla. This approach does not 



require preparation of the system in any subspace since 
all dynamical effects are absorbed by the ancilla. The 
latter is particularly appealing as it avoids the problem 
of imperfect initialization [lo| . The second main result is 
a scheme for universal quantum computation on qubits 
based on this principle. The scheme offers a robust way 
of gate implementation which in comparison to standard 
HQC is less demanding in the preparation stage and in- 
sensitive to certain transitions between energy levels. 

The third main result is a scheme for fault-tolerant 
(FT) HQC based on 2-local interactions. The ultimate 
scalability of any method of computations requires the 
ability for FT implementation that guarantees the exis- 
tence of an accuracy threshold T\\ . The first such HQC 
scheme [l2| requires 3-local interactions that cannot be 
reduced to 2-local using standard techniques without los- 
ing fault tolerance. As 3-local interactions can be very 
difficult to engineer, the question of whether FTHQC can 
be implemented using 2-local Hamiltonians is especially 
important in view of possible experimental implementa- 
tions. We show that our scheme for HQC with noisy 
ancillas is readily compatible with the FT techniques on 
stabilizer codes ll| and can be reduced to 2-local by 
perturbative gadgets [l^] in a FT way. 

Quantum holonomies. — Consider an isodegenerate 
family of Hamiltonians {H\} on an iV-dimensional 
Hilbert space, continuously parameterized by a param- 
eter A in a control-parameter manifold M: H\ — 
Z]n=i Sn(-^)nn(A), wherc {en{X)}n=i ^rc the R differ- 
ent c?„-fold degenerate eigenvalues of Hx, (X]n=i "^n = 
N), and H„(A) are the projectors on the corresponding 
eigenspaces. Let be local coordinates on (1 < M ^ 
AiToM) and A)}^^j^ be an orthonormal basis of the 
71*'' eigenspace of the Hamiltonian at the point A. Con- 
sider a time-dependent Hamiltonian H{t) := Hx{t) ob- 
tained by varying the control parameters along a curve 
A(t) inAi. If the change of A is adiabatic, the transforma- 
tion generated by H{t) is U{t) — Texp(— i dTH(T)) = 
©«^ie''^"(*)C/^Ji), where w„(t) = - dre„(A(r)) are 
dynamical phases, and (t) = Vexp{J^I^^^ ^n), where 
V denotes path-ordering. The adiabatic connections are 



E^^n,^dA^ 



where An^fj, has matrix elements 



(v4„,^)„/3 = {na;X\-^\nl3;X) When the path X{t) 
forms a loop 7(t), 7(0) — ^(T) = Ao, the unitary ma- 
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trix U2 = U^^ (T) = 'Pexp(/^ A„) is called the holonomy 
associated with the loop. The set Hol(A„) = {U2\j e 
La„(A^)}, where L^^^M) = {7 : [0,T] ^ A^|7(0) = 
7(r) = Ao} is the space of all loops based on Aq, is 
a subgroup of U(dn) called the holonomy group. In 
Ref. it was argued that in the generic case, the adia- 
batic connection corresponding to the v}^ energy level of 
the Hamiltonian family is irreducible, i.e., the holonomy 
group Hol(A„) is equal to C/(d„), and therefore adiabatic 
holonomies can be used for universal quantum computa- 
tion in the n"^ eigenspace of the Hamiltonian. 

HQC in subsystems. — A (virtual) subsystem [5] is any 
tensor factor of a subspace of a system's Hilbert space 
H.^ . We will consider decompositions of the form 



(1) 



where logical states are encoded in the subsystems Ti.f', 
and will be interested in universal HQC in {Hf }. Eq. ^ 
describes the most general form of encoding of informa- 
tion and plays a fundamental role in the theory of quan- 
tum error correction 0] where it provides the structure of 
preserved information under open-system dynamics 14 1 . 
(Note that most generally all subsystems Tif can be used 
for encoding and computation simultaneously [l5l|). 

Lemma. Let -ff (0) be a Hamiltonian with at least two 
different eigenvalues. It is possible to implement any 
combination of holonomies in the eigenspaces of i?(0) 
through a suitable adiabatic cyclic change of H{0). 

Comment. It is known that given sufficient control over 
the parameters of a Hamiltonian we can generate holo- 
nomically any unitary operation in a given eigenspace [ij . 
The question we address here is whether it is possible to 
generate an arbitrary combination of holonomies in the 
different eigenspaces. This property may not be obvi- 
ous. For example, in the case of a two-level Hamiltonian, 
the evolution of one eigenspace completely determines 
the evolution of the other one. Since the holonomy in a 
given eigenspace depends entirely on the evolution of that 
eigenspace, this raises the question if the two eigenspaces 
can undergo arbitrary independent holonomies. We now 
show, by construction, that this is possible. 

Proof. It is sufficient to show that it is possible to 
generate a universal set of gates in any given eigenspace 
while at the same time generating the identity operation 
in the rest of the eigenspaces. Without loss of general- 
ity, we will assume that there are only two eigenspaces; 
if there are more, we can always operate within the sub- 
space spanned by two of them by varying only the re- 
striction of the Hamiltonian on that subspace. Then the 
initial Hamiltonian can be written H{0) = eiHi -|- £2^2, 
where Hi^2 are the projectors on the ground and excited 
eigenspaces, and £1 < £2 are their corresponding eigen- 
values. Observe that H(0) is invariant under unitaries of 
the form V = Vi©V2, where Vi 2 are unitaries on the sub- 



spaces with projectors Hi 2, respectively. Let the Hamil- 
tonian vary along a loop H{t), H{0) = H{T), which 
satisfies the adiabatic requirement to some satisfactory 
precision. To this precision, the resulting unitary trans- 
formation can be written U{T) = T exp(-i dtH{t)) = 



'-U2, where Ui and U2 are the holonomies 



resulting in the two eigenspaces, and wi 2 — Jq dtei,2(i) 
are dynamical phases. Observe that the Hamiltonian 
VH{t)V^ , where V — Vi®V2, gives rise to the holonomies 
ViUiVi and V2U2V2 , respectively. This follows from the 
fact that the overall unitary transformation generated by 
VH{t)V^ is equal to VU{t)V'^ where U{t) is the unitary 
generated by H{t). Note that VH{t)V'' is a valid loop 
based on H{0) with the same spectrum as that of H{t). 

Imagine that we want to generate holonomically a gate 
Wi in the ground space of the Hamiltonian while at the 
same time realizing the identity holonomy I2 in the ex- 
cited space. Choose any loop H{t) which gives rise to 
1 

the holonomy in the ground space, where d2 is the 
dimension of the excited space (we know that such a 
loop can be found). Let this loop result in the holon- 
omy W2 in the excited space. The latter can be writ- 
da 

ten T4^2 — X]s*"^li)OI' where {|j)} is an eigenbasis of 

j=i 

W2 and e*"^ , aj G R, are the corresponding eigenval- 
ues. Consider the unitary C2 which cyclicly permutes 
the eigenvectors C2\j) — \j + 1), where we de- 

fine |ci2 + 1) = We can implement the desired 

combination of holonomies in the two eigenspaces as fol- 
lows. First apply H{t). This results in the holonomies 
1 

and W2 in the ground and excited spaces, respec- 
tively. Next, apply {h © C2)H{t){Ii © (72)^ This gener- 
ates the holonomies and C2W2C\ = Y.e^°'^-^\j){j\ 
(we have defined ai-i = a^j). The combined effect of 



these two operations is W-^^ and X] ^ 



i(Qj+aj_i) 



|j)01- We 



holonomies W^^ and C^W2C. 



next apply (/i © Cl)H{t){Ii © C|)t, which generates the 

: |:e"-=|j)(j|. The 

3 d2 

net resuh becomes VF/" and ^ e*("^+"j-i+"^-2)|j)(j|. 

We continue this for a total of 1^2 rounds, which results 

£2 

in the net holonomic transformations W-^^ = Wi and 
gi(ai+a2+...+a<i2) ^ oc /j. This complctcs the proof. 

The proof uses sequences of loops. In the next section, 
we will see that depending on the task it may be possible 
to find constructions based on a single loop. 

Theorem 1. Consider a nontrivial subsystem decom- 
position of the form ([T]). Choose an initial Hamiltonian 
in the form H{0) — ®™ 1 if <S) , where are oper- 
ators on Hf such that all eigenvalues of Hf are different 
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from the eigenvalues of H!^ for i ^ j. In the case when 
m = 1, we impose the additional requirement that 
has at least two different eigenvalues. By varying adia- 
batically this Hamiltonian along suitable loops in a suf- 
ficiently large control manifold, it is possible to generate 
a unitary of the form J7 = 0^ Wf" (g) V^^ , where {Wf"} is 
an arbitrary set of geometric transformations on {Hf}- 

Proof. Denote the eigenvalues of by Wq^, ai = 
l,...,c?i, and the projectors on their corresponding 
eigenspaces Hf. by Ilf. . Then -ff(O) has the spectral 
decomposition iJ(0) = YlT=i I^a,=i ^aj^f ® Hf . , where 
Hf is the projector on Tif-. According to Lemma 1, we 
can implement holonomically any combination of unitary 
transformations in the different eigenspaces of i?(0) up 
to an overall phase. Thus by applying the holonomy 
Wf- ® W^. in each of the eigenspaces Hf ® H^. for = 
l,...,di, where W^. are arbitrary unitaries on Ti^., we 
obtain the net unitary C/ = 0^ W^-^ (g) (0^^ e"''''^W^J = 
0j (g) V^^ , where e*'^"^ are dynamical phases. 

HQC without initialization. — From Theorem 1 one can 
see that in the case when the Hilbert space factors as 
= (8i , it is possible to apply holonomic com- 
putation in subsystem without initializing the state 
of the system in any subspace. In particular, if we are 
given a system Ti.^ in an unknown state, we can append 
to it another ancillary system Ti,^ , also in an unknown 
state, and apply any desired transformation holonomi- 
cally on the first system. Since this approach does not 
require the preparation of pure ancillary states, it can 
be advantageous in implementations where the latter is 
difficult, such as nuclear magnetic resonance (NMR) [l^ . 

We now present an explicit scheme for universal com- 
putation on qubits based on this principle. The scheme 
uses a single ancillary gauge qubit. We will show how to 
implement a universal set of one- and two-qubit gates. 
We will build the necessary loops by interpolations be- 
tween points in the space of Hamiltonians of the form 
H{t) = f{t)H{0) + 9{t)H{T), where /(O) = g{T) = 1, 
/(T) ^ (7(0) = 0. This Hamiltonian interpolates between 
i?(0) and H{T) during a time interval T. The interpolat- 
ing Hamiltonians that we will be using (to be described 
below) have two energy levels of equal degeneracy, and 
their energy gaps are non-zero unless the entire Hamilto- 
nian vanishes, i.e., the interpolations can be realized adi- 
abatically for sufficiently long time T and smooth choice 
of f{t) and g{t). For an adiabatic interpolation of the 
above type we will use the short notation H{0) H{T). 

Let us label the two qubits on which we will be ap- 
plying the gates by 1 and 2, and the gauge qubit by 3. 
In order to apply a single-qubit gate, say, on qubit 2, we 
will use the starting Hamiltonian H{0) = I2® X^,. (Here 
by Xj, Kj, Zi and li we denote the Pauli matrices and 
the Identity acting on the i^^ qubit). We first apply 



17|. 



(2) 



This results in a geometric transformation and a dy- 
namical phase in each eigenspace. Let us denote the 
purely geometric part of the resulting unitary by C/2,3 
(the exact form of C/2,3 is not important since we will 
undo it later). At this point we can apply a unitary 
whose geometric part is equal to an arbitrary gate on 
qubit 2 according to a method described in Ref. 
For example, the interpolation Z2 ® Z^ X2 
gives rise to the geometric operation R2Z2, where R2 
is the Hadamard gate on qubit 2. The interpolation 
Z2(^Z3 -(cosf X2-hsin|y2)«'Z3 -> -Z2<»Z3 results 
in the geometric operation T2Z2X2 , where T2 denotes 
the 7r/8 gate on qubit 2. These two gates are a universal 
set of single-qubit gates. Let G2 be the gate from the 
above set which we want to implement. After the corre- 
sponding interpolation, the net geometric transformation 
becomes G2C/2,3 and the Hamiltonian is transformed to 
G2Z2G\^Z3. We can now "undo" the unitary C/2,3 by ap- 
plying the interpolation G2.Z2G2<8'Z3 12^X3. The lat- 
ter is the inverse of Eq. ^ up to the single-qubit unitary 
transformation G2, i.e., it results in the transformation 
G2UI3GI. Thus the net result is G2C/2%GJG2 C/2,3 = G2, 
which is the desired unitary on qubit 2. Note that the 
relative dynamic phase between the ground and excited 
spaces, which accumulates during the procedure, at the 
end is equivalent to a transformation on qubit 3. 

For universal computation, we also need a nontrivial 
two-qubit gate. We can start again by the interpolation 
(PI which results in the geometric transformation C/2,3. 
At this point we can apply, for example, the interpolation 
Ii(g) Z2(E) Z3 -> /i (g) I2 ® ^3 Zi® Z2® Z-i, which 
results in the gate S'|Afi,2, where 5*1 = and /V1.2 is the 
"controlled not" gate with qubit 1 the control, and qubit 
2 the target To "undo" the operation C/2,3, we apply 
the transformation Z^® Z2® Z^ ^ 1\®I2® ATa, which is 
the inverse of ^ up to the transformation S'Ja^i,2. The 
net result is 5|Ari,2C/|,3/Vi,25i5|/Vi,2C/2,3 = S\Xx^2- 

We note that unlike the standard holonomic approach, 
here each eigenspace of the Hamiltonian undergoes the 
same geometric operation, which supplies the scheme 
with additional robustness. The scheme is insensitive to 
those transitions between the two energy levels that are 
equivalent to local operations on the transformed gauge 
qubit. The scheme uses 2- and 3-local Hamiltonians. 

FTHQC with 2-local Hamiltonians. — The theory of 
fault tolerance ll| guarantees that, if errors during the 
implementation of a given gate are sufficiently uncorre- 
lated and improbable, an arbitrarily long computation 
can be implemented reliably with a modest resource over- 
head. The first proposal for FTHQC uses the encod- 
ing present in a stabilizer code and Hamiltonians that are 
elements of the instantaneous stabilizer or gauge group of 
the code. These Hamiltonians couple qubits in the same 
block, but errors do not propagate as each eigenspace is 
subject to the same transversal operation. That scheme 
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requires 3-local Hamiltonians. Even though every Hamil- 
tonian can be simulated by a 2-local one via the so called 
perturbative gadgets [l3| , the locality of that scheme can- 
not be reduced by a direct application of these tech- 
niques. Since the simulated Hamiltonian couples qubits 
in the same code block, an error on one of the gadget an- 
cillas [3] can spread to multiple qubits within a block. 

The scheme from the previous section suggests an al- 
ternative approach to FTHQC. Any FT protocol on qubit 
stabilizer codes can be decomposed into transversal one- 
and two-qubit gates (these are gates that couple only 
the corresponding qubits from different blocks) . In addi- 
tion, one requires the preparation of a special ancillary 
state such as (|00...0) -I- |ll...l))/-\/2, which is done non- 
transversally. Transversality guarantees that a single er- 
ror during an encoded operation results in at most one 
error per block of the code. Our scheme for holonomic 
one- and two-qubit gates is readily compatible with this 
approach: we can apply the same operations as in a stan- 
dard FT protocol [Ul by coupling every qubit or pair of 
qubits in the code to an "external" gauge qubit as de- 
scribed in the previous section. Obviously, a single er- 
ror during the implementation of a transversal operation 
cannot propagate to multiple qubits in a block because 
the latter do not interact. In contrast to the previous 
approach which can be understood as performing HQC 
inside the subsystem containing the protected informa- 
tion, this scheme performs HQC in the entire system and 
does not require Hamiltonians that depend on the code. 
Furthermore, here each 3-local Hamiltonian can be re- 
duced to 2-local via perturbative gadgets as the gadget 
ancillas would couple to at most one qubit inside a block. 
A complete fault-tolerance analysis is beyond the scope 
of this paper, but we note that the use of extra qubits 
increases the chance for an error during a single gate. In 
addition, the 3-qubit gadget decreases the minimum gap 
of the original Hamiltonian by a factor where e is 

the perturbation parameter in which the approximation 
to order 0(£^) is carried out [l^. Thus for maintaining 
a given precision, the time for implementing a two-qubit 
gate would have to increase by a similar factor which de- 
creases the allowed rate for environment noise. A way 
around this could be to look for non-perturbative imple- 
mentations. Since the no-initialization property is not 
crucial for fault tolerance per se, implementations with 
an ancilla in a known state can also be considered. 

Conclusion. — In summary, we have introduced a gen- 
eral framework for HQC in subsystems, showing that it 
is possible to realize simultaneously independant HQC 
in the subsystems Tif in any nontrivial decomposition 
of the form ([1]). As an apphcation, we proposed a ro- 
bust scheme for applying purely geometric gates to un- 
encoded qubits by the use of a noisy ancillary qubit. 
We used this approach to show that 2-qubit Hamiltoni- 
ans are universal for FTHQC. We hope that our results 
will open new avenues for quantum information process- 



ing implementations that combine the robustness of the 
holonomic control with the most general form of encod- 
ing. An interesting future direction would be to extend 
the present results to the theory of geometric phases 
based on dynamical invariants 18|, which encompasses 
non-adiabatic, mixed-state and open-path holonomies. 
Since geometric phases have wide applications, associ- 
ating holonomies with subsystems could find use beyond 
the field of quantum computing as well. 
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Note. — Recently, an independent scheme for (open- 
path) FTHQC with 2-local Hamiltonians via perturba- 
tive gadgets was proposed by D. Bacon and S. T. Flam- 
mia in Ref. [19i |. 
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